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Abstract. We present a long (150 ks elapsed time) X-ray ob- 
servation of the dM3.5e star EV Lac, performed with the ASCA 
observatory in July 1998, during which an exceptionally in- 
tense flaring event (lasting approximately 12 ks) was observed: 
at the flare's peak, the X-ray count rate in the ASCA GIS de- 
tectors was ~ 300 times the quiescent value. The physical pa- 
rameters of the flaring region have been derived by analyzing 
the decay, using both a "classic" quasi-static approach and an 
approach based on hydrodynamic simulations of decaying flar- 
ing loops. Notwithstanding the large peak X-ray luminosity, 
this second method shows that the flare's decay is compatible 
with its being produced in a relatively compact region of semi- 
length L ~ 1.3 X 10^° cm (~ 0.5 R^), large but not excep- 
tional even by solar standards. The flare decays is fast (with a 
measured e-folding time for the light curve of < 2 ks), but nev- 
ertheless the hydrodynamic -based analysis shows strong evi- 
dence for sustained heating, with the shape of the light curve 
dominated by the time evolution of the heating rather than by 
the natural cooling of the flaring plasma. As a consequence, the 
quasi-static method yields a much larger estimate of the loop's 
length (L ~ 2 i?*). The event shows (similarly to some other 
well-studied large stellar flares) a light curve characterized by 
two separate decay time constants (with the initial decay be- 
ing faster) and a significant enhancement in the plasma metal 
abundance at the peak of the flare. The energetics of the event 
are exceptional, with the peak X-ray luminosity of the event 
reaching up to ~ 25% of the bolometric luminosity of the star, 
making this one of the largest X-ray flare (in relative terms) 
observed to date on a main-sequence star. 

Key words: Stars: individual: EV Lac - Stars: late-type - Stars: 
activity - Stars: coronae; X-rays: stars 



1. Introduction 

One of the basic standing problems of stellar coronal astron- 
omy is the determination of the spatial structuring of the coro- 
nal plasma. While for the study of stellar interior structure the 
first-order approximation of spherical symmetry is a good start- 



ing point, stellar coronae are, as shown by the large body of ex- 
tant imaging observations of the solar corona, far from spheri- 
cally symmetric. The solar corona shows a high degree of spa- 
tial structuring: most of the X-ray luminous plasma is con- 
fined in coronal loops preferentiaUy located at mid-latitudes, 
with an average position which tracks the migration of sunspots 
through the solar cycle. The lack of spatial information consti- 
tutes a strong limitation for the study of stellar coronae: low- 
resolution coronal X-ray spectra are insensitive to the plasma's 
density, so that non-dispersive, CCD- or proportional counter- 
based spectral observations do not allow to distinguish between 
a large diffuse corona and a compact, structured, high-pressure 
one. It is thus not possible, if the solar analogy is postulated, 
to determine how the solar corona scales toward higher levels 
of activity, i.e. if largely through a spatial filling of the avail- 
able volume with coronal loops (yielding a relatively symmet- 
ric corona) or if through the filling of a relatively small number 
of coronal structures with significantly higher density plasma. 

Thus far the main tools to study the spatial distributions 
of the coronal plasma have been eclipse experiments and the 
study of flares. While the observation of rotational modulation 
should also in principle allow to derive t he spatial distrib ution 
of the emitting plasma, as discussed by Schmitt (1998)| con- 
vincing examples of rotationally modulated X-ray emission are 
rare. Under a given set of assumptions, the study of the decay 
phase of a flare can yield information about the size of the flar- 
ing structure; different methods for this type of analysis have 
been developed and applied in the past to several observations 
of stellar flares. The widely applied quasi-static approach (see 
below) almost invariably results, when applied to intense stellar 
flares, in long, tenuous coronal loops extending out to several 
stellar radii. The stronger flares yield in general larger sizes. 
Detailed hydrodynamic modeling of flaring loops has provided 
useful insight on stellar X-ray flares ( ^eale et al. 1988 ); more 
recently, diagnostic tools have been developed for the deriva- 
tion of the size of stellar flaring loops and of the heating evo- 
lution dReale et al. 19931; iReale et al. 19971; Reale & Micela 
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1998). In the solar case, in addition to the "compact" flares, 
in which the plasma appears to be confined to a single loop 
whose geometry does not significantly evolve during the event 
(an assumption shared by the quasi-static method and by the 
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hydrodynamic modeling quoted above), a second class of flares 
is usually recognized, i.e. the "two-ribbon" events, in which 
an entire arcade of loops first opens and then closes back; the 
footpoints of the growing system of loops are anchored in Ha- 
bright "ribbons". These flares are generally characterized by 
a slower rise and decay, and a larger energy release. Compact 
flares have often been considered to be due to "impulsive" heat- 
ing events, while the longer decay times of two-ribbon events 
have been considered as a sign of sustained heating. However, 
also in the case of compact flar es sustained heati ng has been 
shown to be frequently present ( |Reale et al. 1997 ), so that the 
distinction may indeed be less clear than often thought. 

Fitting of X-ray spectra with physical models of static loops 



flares come in a large variety of sizes for both their time scales 



and their energetics. Pallavicini et al. (1990) did not attempt to 



(Maggio & Peres 1997; Sciortino et al. 1999) can also yield 



the surface filling factor of the plasma as one of the fit pa- 
rameters. However, for loops smaller than the pressure scale 
height this method only gives an upper limit to the filling 
factor, which needs to be further constrained, for example, 
with estimates of the pl asma density based on EUV line ratios 
( Maggio & Peres 1997 ). With few exceptions[|, eclipse studies 
of the quiescent emission of eclipsing binaries have thus far 
failed to yield strong c onstraints on spatial structuring of the 
plasma ( Schmitt 1998 ). In part this is due to the inversion of 
the weak observed modulation being mathematically an intrin- 
sically ill-posed problem, where few compact structures can 
mimic the emission from a more diffuse medium (see discus- 



sion in Schmitt 1998). Recently, the observation of the total 
eclipse of a large flare on Algol (Schmitt & Favata 1999) has for 
the first time yielded a strong geometrical constraint on the size 
of a flaring structure. The geometrical loop size is significantly 
smaller than the size derived from the analysis of the flare's 
decay ( favata & Schmitt 1999| ) using the quasi-static method, 
showing how such approach can over-estimate the actual loop 
size. The characteristics of the Algol flare are such that the hy- 
drodynamic decay, sustained heating framework which we also 
use here yields a large range of allowed loop semi-lengths, with 
the lower end of the range marginally compatible with the ge- 
ometrically derived size. 

The presence of intense X-ray flares on flare stars^ was 
well established with Einstein observations - with some by 
now classic observations such as the one relative to a flare on 
Proxima Cen ( ^aisch 1983 ). However, Einstein observations 
were usually relatively short (few ks) thus imposing a bias on 
the type of events which could be detected (Ambruster et al. 
1987). In particular, the short total integration times reduced 
the chance of detecting rare event types. The EXOSAT obser- 
vatory featured long, uninterrupted observations, and thus al- 
lowed to collect a more unbiased view of flares on low-mass 



stars, resulting in a database ( Pallavicini et al. 1990 ) of about 
300 hr of flare stars observations, from which it is apparent that 



' Notably the observation of a CrB, Schmitt & Kurster (1993j , in 
which the observed total eclipse constrains the corona on the G5V star 
to have a scale height much less than a solar radius. 

^ We use the term to mean "M-type dwarfs, either single or mem- 
bers of a multiple system, which show frequent sudden enhancements 
of their X-ray, UV and optical luminosity". 



derive the spatial scales of the observed flares, although they 
tentatively divided them into two classes reminiscent of the 
solar compact and two-ribbon flares. The next generation of 
of flare observations came with the ROSAT All-Sky Survey 
(RASS), which, thanks to its sky scanning strategy, allowed 
to search for rare, long-lasting events. Some flares of previ- 
ously undetectedmagnitude and duration are indeed present 
in it ( Schmitt 1994 ); in particular, EV Lac showed a long X- 
ray flare lasting approximately a day, superimposed on a much 
shorter but more intense event. 

We have observed EV Lac for two days with the ASCA 
observatory, detecting the most intense X-ray flare thus far ob- 
served on a main-sequence star, with a 300-fold peak increase 
of the X-ray count rate. This paper presents a detailed analysis 
of the flaring event, and it is organized as follows: the ASCA 
observations and their reduction, together with the determina- 
tion of the spectral parameters for the flare are discussed in 
Sect. H the parameters of the flaring region are determined (us- 
ing both the quasi-static formalism of van den Oord & Mewe 
1989 a nd the hydrodyna mic decay, sustained heating frame- 
work of Reale et al. 1997 ) in Sect. ||, with a discussion (Sect. Q) 
closing the paper. 

2. Observations and data reduction 



EV Lac was observed by ASCA continuously for ~ 150 ks 
from 13 July 1998 06:10 UT to 15 July 1998 01:40 UT, as an 
AO-6 guest investigator target (P.I. F. Favata). The data were 
analyzed using the FTOOLS 4.1 software, extracting both the 
spectra and the light curves with the XSELECT package and 
performing the spectral analysis with the XSPEC 10.0 pack - 
age. The MEKAL plasma emission model ( [VIewe et al. 1995 ) 
was adopted for the spectral analysis. The peak count-rate of 
the flare (~ 100 cts s^^ in the SIS) exceeds both the teleme- 
try limit (~ 40 cts s~^) and the 1% pile-up limit throughout 
the whole point-spread function (~ 60 cts s^^), thus prevent- 
ing a reliable spectroscopic analysis. We therefore only used 
the GIS data in the following. Source photons have been ex- 
tracted, for both GIS-2 and GIS-3 detectors, from a circular 
region 18.5 arcmin in diameter (74 pixels, somewhat larger 
than the suggested extraction radius for GIS data; given the 
strength of the source this ensures that as many as possible 
source photons are collected) centered on the source position, 
while background photons have been extracted from a circular 
region identical in size to the source region but symmetrically 
placed with respect to the optical axis of the X-ray telescope. 
For point sources such strategy allows for the background to be 
extracted from the same observation (and thus with the same 
screening criteria) while ensuring that the effect of telescope vi- 
gnetting on cosmic background photons is properly accounted 
for Given the high intensity of the source emission during the 
flare, the background is however effectively negligible. The 
GIS-3 background-subtracted light curve for EV Lac for the 
complete duration of the ASCA observation is shown in Fig. |l]. 
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Fig. 1. The background-subtracted (background count rate ~ 0.05 cts s ^) GIS-3 light curve of EV Lac for the entire ASCA 
observation, binned in 180 s intervals. 
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Fig. 2. The background-subtracted (background count rate ~ 0.05 cts s^^) GIS-3 light curve of the large EV Lac flare discussed in 
the present paper, binned in 60 s intervals. The extent of the time intervals in which the flare has been divided for the determination 
of the flare's spectral parameters is shown. No dead time correction has been applied to the count rates. 



The light curve shows evidence for variability on many 
time scales, and at least three individual flares can be identified: 
the exceptional event at ~ 51 ks from the beginning of the ob- 
servation, and two minor (however still rather sizable) flares at 
~ 10 and ~ 75 ks. The light curves of both minor events show 
a clear exponential decay, but their peak is not covered by the 
observations. To derive the temporal evolution of the tempera- 
ture and emission measure of the large flare we have subdivided 
it in 9 time intervals, shown together with the GIS-3 light curve 



of the event in Fig. ^. Individual GIS-2 and GIS-3 spectra have 
been extracted for each of these intervals and merged using the 



procedure described in the ASCA ABC Guide (1997). The ex 



posure time of each individual spectrum has been corrected for 
the dead-time of the GIS (which, at these high count rates is 
rather significant, with values up to 1.2 - note that the both 
light curves from Figs. |l] and ^ are not corrected for detector 
dead-time). The quiescent spectrum has been taken from the 
interval immediately preceding the flare (interval 1 in Fig. |[ 
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covering ~ 5 ks). A two-temperature model has been fit to the 
spectrum extracted from interval 1, with the resulting spectral 
parameters shown in Table |l]. 

Table 1. The spectral parameters derived for the quiescent 
emission of EV Lac from the analysis of the GlS-2 and GlS-3 
spectra accumulated during the time interval 1, using a two- 
temperature MEKAL spectral model. The fit has 72 degrees of 
freedom, and the corresponding probability level is 99.8%. The 
quiescent X-ray luminosity corresponding to the above spectral 
parameters is Lx = 3 x 10^^ erg s^^. 
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The sequence of GIS flare spectra is shown in Fig. g|. 
The left panel shows the spectra collected during the rising 
phase of the flare, up to the peak in the light-curve (inter- 
val 4), while the decay-phase spectra are plotted in the right 
panel. Individual flare spectra from time intervals 2 to 9 have 
been fit with a single-temperature MEKAL model; given the 
lack of soft response in the GIS no absorbing column density 
was included. The quiescent emission was modeled by adding 
a frozen-parameter two-temperature MEKAL model to the fit 
(with the parameters as in Table |l]); the results of the spec- 
tral fits are shown in Table ^. During intervals 5 and 6 the 
single-temperature fit does not yield a satisfactory reduced 
(see Table ||), due to the large residuals present in the region 
around 1 keV, with the observed spectra showing some addi- 
tional emissivity with respect to the models. Similar effects are 
also seen in the in the ASCA SIS Capella spectrum (Brickhouse 
1998) and in the flaring spectra of Algol as seen by SAX (Fa- 
vata & Schmitt 1999), suggesting that current plasma emission 
codes (as the MEKAL one used here) under-predict the observed 
spectrum around 1 keV, likely due to a large number of high 
quantum-nu mber (n > 5) Fe L lines from Fe XVII, Fe XVIII 
and Fe XIX ( Brickhouse 1998 ), and thus the formally unaccept- 
able resulting from the fit does not necessarily imply that 
the one-temperature model is not correctly describing the ob- 
served spectrum. Indeed, we have verified that adding further 
temperature components does not significantly improve the fit 
for intervals 5 and 6. The time evolution of the flare's spec- 
tral parameters (temperature, emission measure, plasma metal 
abundance) is shown in Fig. ^, together with the flare's GIS-3 
light-curve binned in the same time intervals. 

3. Flare analysis 

We have analyzed the present flare using two different frame- 
works (quasi-static cooling and hydrodynamic decay, sustained 
heating), which both make the assumption that the flaring 
plasma is confined in a closed loop structure, whose geometry 
is not significantly evolving during the event. Although direct 



support to this assumption is of course missing, the relatively 
short duration of the event allows an analogy with solar com- 
pact flares. In any case the second method provides reliable 
scale sizes of the flaring structure even in the presence of some 
readjustment of the magnetic field, the crucial assumption been 
plasma confinement. 

3.1. The quasi-static cooling framework 

Many stellar X-ray flares observed to date have been studied 
using the quasi-static formalism first discussed in detail by 
van den Oord & Mewe (1989) It is thus of interest to ana- 



lyze the present event with the quasi-static approach, to allow 
a homogeneous comparison with literature data, even if, as dis- 
cussed by favata & Schmitt (1999) , this method can signifi- 
cantly over-estimate the size of the flaring loops (see also Reale 
et al. 1997). 

According to van den Oord & Mewe (1989)| , for the decay 
to be quasi-static (i.e. to happen through a sequence of states 
each of which can be described by the scaling laws applicable 
to stationary coronal loops) the ratio between the characteris- 
tics times for radiative and conductive cooling must be constant 
during the flare decay (although its absolute value needs not be 
known). This ratio is parameterized by the quantity 



A* = 



C X 



rpl3/i 

EM ' 



(1) 



The normalization C depends on the details of the loop's 
geometry, and is not relevant here. The evolution of fi is plotted 
in Fig. U within the error bars /i is constant during the whole 
decay, so that the conditions for the applicability of the quasi- 
static decay framework are in this case met. 

The scaling laws discussed by Stern et al. (1992^ -linking 



the loop semi-length L and the plasma density n with the peak 
flare temperature T and the effective decay time t - yield re- 
sults very similar to the ones obtained through a full fit to the 



equations of van den Oord & Mewe (1989), so that we will 
limit ourselves to their application. Specifically, L cx tT'^/^ 
and n cx r^^T^^^, valid for temperature regimes in which the 
plasma emissivity scales as ^/qT'"^ with 7 ~ —0.25. In prac- 
tice this is valid for T ^ 20 MK, i.e. for most of the decay of 
the flare discussed here. Scaling the values determined for the 
EV Lac flare from the parameters derived by v an den Oord & 
Mewe (1989) for the EXOSAT flare observed on Algol, the de- 
rived loop semi-length is L ~ 5 x 10^" cm (~ 2 i?*), and the 
plasma density n : 



3.2. The hydrodynamic decay, sustained heating framework 

A different approach to the study of a flare's decay phase, with 
the same aim of determining the physical parameters of the flar- 



ing region, has been developed by Reale et al. (1997). It has 
been recognized from the modeling of solar X-ray flares that 
the slope of the locus of the flare decay in the log n-log T plane 
is a powerful diagnosti c of the presence of a dditional heating 
during the decay itself ( Sylwester et al. 1993); by making use 
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Fig. 3. The time sequence of merged (GIS-2 + GIS-3) spectra for the EV Lac flare. The left panel shows the sequence of spectra 
during the rise phase of the flare, up to its peak, while the right panel shows the sequence of decay spectra. The numeric label 
at the left of each spectrum indicates the time interval during which the spectrum was collected, as numbered in Fig. |[ Each 
spectrum is overplotted with the corresponding best-fit model from which the spectral parameters have been derived. The spectra 
have been rebinned to a minimum signal-to-noise ratio per bin of 2. 



Table 2. The spectral parameters T, Z and EM derived for the individual phases of the EV Lac flare from the analysis of the GIS 
spectra accumulated during the time intervals 2 to 9. The spectra have been analyzed with a single-temperature MEKAL model 
(plus a frozen-parameter two-temperature model to account for the quiescent emission). The bounds of the confidence intervals 
at the 90% levels (Ax^ = 3.5) are also reported for each parameter. The time interval i to which each set of parameters applies 
is shown in Fig. ^. 
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of extensive hydrodynamic modeli ng of decaying flari ng loops, 
with different heating time scales, |^eale et aT (1997)1 have de- 
rived an empirical relationship between the light curve decay 
time (in units of Tth, the loop thermodynamic decay time, Se- 
rio et al. 1991) and the slope C in the log n-log T diagram. This 
allows to derive the length of the flaring loop length as a func- 
tion of observable quantities, i.e. the decay time of the light 
curve, the flare maximum temperature and the slope of the de- 
cay in the log n-log T diagram (the square root of the emission 
measure of the flaring plasma is actually used as a proxy to the 
density). Since the characteristics of the observed decay depend 
on the specific instrument response, the parameters of the ac- 



tual formulas used have to be caUbrated for the telescope used 
to observe the flare. 



The method reported in |leale et al. (1997) was tested on 
a sample of solar flares observed with Yohkoh-SXT, for which 
both images (from which the length of the flaring loop was 
measured) and spectral parameters (from the temperature and 
emission measure diagnostic derived from Yohkoh filter ratios) 
were available, and has been shown to provide reliable results 
for most of the studied events. A first application of the method 
to stellar flares observe d with ROSAT PSPC is described by 
Reale&Micela(1998). 
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Fig. 4. The temporal evolution of the spectral parameters of the flare, i.e. the plasma metal abundance (top left), emission measure 
(top right) and temperature (bottom left). Also shown (bottom right), the light-curve of the flare as seen by the GlS-3 detector, 
integrated in the same temporal intervals as the ones used to extract the flare's spectra. The best-fit exponential decay to the light 
curve is also plotted. The dashed horizontal line in the abundance plot represents the value of the metal abundance determined 
for the quiescent, pre-flare spectrum. 



For the present study the method has been recalibrated for 
stellar flares observed with ASCA GIS. The thermodynamic 
decay time rth of a closed coronal loop w ith semi-length L 
and maximum temperature T^ax is given by Serio et al. (1991) 
as 



Tth 



aL 

vTma 



(2) 



where a = 3.7 x 10"^ cm~^ K^/^. By means of a grid of 
hydrostatic loop models (see |^eale & Micela 1998 ) we have 
determined an empirical relationship linking the loop maxi- 
mum temperature Tmax^ typically found at the loop apex (e.g. 



Rosner et al. 1978 ) to the maximum temperature Tobs deter- 
mined from the GIS spectrum: 



= 0.085 x Tib^^^ 



(3) 



Following the same procedure as in Reale et al. (1997) and 



Reale & Micela (1998) (using extensive hydrodynamic mod- 



eling of decaying flaring loops) we have determined the ratio 
between tlc (the observed e-folding time of the flare's light 
curve determined by fitting the light curve from the peak of 
the flare down to the 10% of peak level) and rth as a function 
of the slope C in the log V EAd -log T diagram to be, for the 
ASCA GIS 



— = F(C) 

Tth 



(4) 
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Fig. 5. The temporal evolution of the ji parameter (the ratio be- 
tween the radiative and conductive cooling time for the loop) 
during the flare. The dotted line joins the points corresponding 
to the different time intervals reported in Fig. |[ 



where Ca = 10.9, Ca — 0.56 and qa — 0.6. The formula for the 
loop semi-length L is therefore: 



TLCvTma 



0.38 < C < 1-7 



(5) 



where the second part of the relationship gives the range of ^ 
values allowed according to the modeling. The uncertainty on 
L is the sum of the propagation of the errors on the observed 
parameters tlc and with the standard deviation of the differ- 
ence between the true and the derived loop lengths. The latter 
amounts to ~ 15%. Equation ^ has been tuned on exponentially 
decaying light curves; however it has been shown to provide 
reliable results also on solar flares with more complex decay 
trends, e.g. a double exponential decay (as for the flare stud- 
ied here), provided that the whole decay is fitted with a single 
exponential (F. Reale, private communication). 

The evolution of the EV Lac flare in the log V EM-\og T 
plane is shown in Fig. |[ together with a least-square fit to 
the decay phase. The resulting best-fit slope for the decay- 
ing phase computed for time intervals from 4 to 9 inclusive 
is ^ = 0.56 ± 0.04. Application of Eq. ^ yields a ratio between 
the observed cooling time scale tlc and the thermodynamic 
cooling time scale for the flaring loop Tth of F{C) ~ 4.6. Such 
a large value implies that the observed decay is driven by the 
time-evolution of the heating process and not by the free cool- 
ing of the loop. Also, the actual loop length will be significantly 
smaller than it would be estimated assuming that the spectral 
parameters reflect free cooling of the decaying loop. The ac- 
tual value of Tlc has been determined by fitting the GIS light 
curve, binned in the same time intervals used for extracting the 




Fig. 6. The evolution of the flare's decay in the log V EM- 
logT plane. The points plotted represent the flare's evolution 
from time interval 2 to time interval 9 inclusive. The dotted line 
joins the points corresponding to successive intervals. The de- 
cay phase begins with the third interval, and it closely follows 
a straight line. The continuous line is a least-square fit to the 
decay phase, with a best-fit slope = 0.56 ± 0.04, while the 
numeric labels indicate the time interval to which each point 
refers. 



individual flare spectra (as plotted in Fig. Q), considering the in- 
tervals from 4 to 9 inclusive. In this case tlc — 1-80 ± 0.15 ks, 
and therefore Tth — 400 s. The intrinsic flare peak temper- 
ature is, applying Eq. || to the observed maximum tempera- 



ture, Tmax - 150 MK. From Eq. 
L = (1.3 ±0.3) X 10^° cm, i.e. L 



the loop semi-length is 
0.5 i?*. This loop length 



is much smaller than the pressure scale height^ of the flaring 
plasma on EV Lac and also significantly smaller (by a factor of 
4) than the one derived through the quasi-static formalism. 

A simple consistency check can be obtained by comparing 
the pressure obtained by assuming that the flaring loop is not, at 
maximum, far f rom a steady-state c ondition (thus applying the 
scaling laws of Rosner et al. 1978 ) with the pressure implied 
by the derived values of L and T for a plausible geometry. In 
practice the geometry is parameterized by the ratio /3 between 
the radius of the loop and its lenght. The pressure is then 



/ EM 



(6) 



If we assume (3 ~ 0.1-0.3 (a typical range for solar coro- 
nal loops) the loop volume is ~ 1.4-13 x 10^^ cm'^, and the 

defined as H — 2kT/fig, where T is the plasma temperature in 
the loop, /i is the molecular weight and g is the surface gravity of the 
star. In this case iif ~ 8 x 10^^ cm 
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resulting plasma density at the peak of the flare is n ~ 2- 
0.2 X 10^^ cm^"^. The corresponding pressure is p,„ax — 8- 
0.9 X 10'' dyne cm^^. Using the scaling laws of Rosner et al. 
(1978) applicable to steady-state loops, 



Tlnax = 1.4 X 103(poi)'/^ 



(7) 



where pq is the pressure at the base of the loop, one obtains 
Po ~ 10^ dyne cm^^, slightly larger than Pmax for /3 = 0.1. 
This implies that the plasma evaporated from the chromosphere 
has not (/3 = 0.3) or nearly {[3 = 0.1) filled the flaring loop up 
to the hydrostatic equilibrium conditions at flare maximum. 

3.3. Energetics 

We have computed, for each of the 8 time intervals in which 
the flare has been subdivided, the X-ray luminosity in the 0.1- 
10. keV band. For this purpose the spectrum has been extrap- 
olated outside of the formal spectral range covered by the GIS 
detectors; this is at most likely to underestimate the true lumi- 
nosity in the extended band, as it may miss any softer com- 
ponent present in the spectrum and invisible to the GIS. The 
time-evolution of the flare X-ray luminosity is shown in Fig. ^ 
in which the data are plotted both in absolute units and in units 
of the star's bolometric luminosity. 



J I , , I 1 I , , , , I , , , , I , , L 

4.5xl0« 5x10* 5.5x10* 6x10* 

Seconds since UTC 30/08/97 03:16:23 

Fig. 7. The temporal evolution of the flare's X-ray luminosity 
calculated in the 0. 1 to 10. keV band starting from the quiescent 
state. The left axis gives the luminosity in erg s^^, while the 
right axis shows the same value in units of the star's bolometric 
luminosity (5.25 x lO'^^ erg s^^). 



During interval 4 (at the peak of the light curve) the X-ray 
luminosity of the flare is about one quarter of the photospheric 
(bolometric) luminosity of the star. Soft X-ray radiation is only 



one of the energy loss terms for the flaring plasma, with kinetic 
energy, conduction to the chromosphere and white light, UV 
and XUV flaring emission also significantly contributing to the 
energy budget. In the solar case, detailed analyses of flares of 



different types (Wu et al. 1986) show that at the peak of the 



event soft X-ray radiation only accounts for 10-20% of the total 
energy budget; similarly, |ioudebine et al. (1993) analyzed a 
large optical flare on the dMe star AD Leo, concluding that the 
kinetic energy of plasma's motions during the event is likely to 
be at least as large as the radiated energy during the flare. 

A detailed assessment of the energy balance of the present 
flare is not possible, given the lack of multi-wavelength cover- 
age and of velocity information which could help assessing the 
plasma kinetic energy. The total energy radiated in the X-rays 
is E ~ 1.5 X lO^'' erg (obtained with simple trapezoidal inte- 
gration of the instantaneous X-ray luminosity), over ~ 10 ks, 
equivalent to ~ 300 s o f the star's bolometri c energy output. 
From the scaling laws of Rosner et al. (1978) we can estimate 



the heating rate per unit volume at the peak of the flare, as- 
sumed uniform along the loop, as 



dH 
dVdt 



10^ p'^/^ L-^/^ ~ 240 erg cm'^ s" 



(8) 



The total heating rate at the flare maximum is therefore 



dH 
~dt 



dH 
dVdt 



X V 



3.3 X 10^^ ergs-i 



(9) 



a factor of ~ 3 higher than the flare maximum X-ray luminos- 
ity (see Fig. compatible with X-ray radiation being only one 
of the energy loss terms during the flare. If we assume that the 
heating is constant for the initial rising phase, which lasts for 
tf ~ 300 s, and then decays exponentially, with an e-folding 
time th — 4.6Tth — 1800 s (i.e. similar to tlc), the total en- 
ergy released during the flare is 



Ht, 



dH 
~dt 



X {tr 



7 X 10^** 



erg 



(10) 



approximately five times the energy radiated in X-rays. Energy 
losses by thermal conduction are indeed expected to be large at 
such high temperatures. 

3.4. Metal abundance 

The metal abundance of the flaring plasma shows a significant 
evolution, rising, at the peak of the flare, to a value ~ 2 times 
higher than the abundance measured for the quiescent emis- 
sion, and decaying back to the quiescent value during the ter- 
minal phase of the flare. Evolution of the metal abundance has 
been observed also in other flares - for example it is evident in 



the Algol flare observed with SAX, Favata & Schmitt (1999) 



and the behavior observed here appears to follow the same gen- 
eral pattern of abundance enhancement going in parallel with 
the flare's light curve. The simplest explanation suggested for 
this has been to assume that a fractionation mechanism is at 
work in the quiescent corona, and that the evaporation of pris- 
tine photospheric material during the flare temporarily sets the 
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coronal chemical equilibrium off balance. Unfortunately, few 
reliable abundance determinations are available for M dwarfs; 
in particular, no state-of-the-art photospheric abundance anal- 
ysis of EV Lac is known to us, and thus it is not possible to 
assess whether the quiescent coronal abundance is indeed de- 
pleted with respect to the photospheric value. 

4. Discussion 

The most remarkable characteristic of the EV Lac flare dis- 
cussed in the present paper is certainly its very large X-ray lu- 
minosity: at flare peak, for a few minutes, the flare nearly out- 
shines the star's photosphere. A detailed analysis of the flare's 
decay, however, shows that this is an interesting event also on 
other accounts. Typical loop lengths derived for strong flares on 
active stars, mostly using the quasi-static cooling mechanism, 
are large, comparable or often greater than the stellar radius. 
The picture which has emerged from most of the literature is 
thus one of long, tenuous plasma structures, with the attendant 
challenges of sufficiently strong magnetic fields far away from 
the stellar surface. In the present case, instead, the loop semi- 
length deri ved from the analys is of the flare decay using the 
method of Reale et al. (1997) is relatively compact, at about 
0.5 stellar radii (implying a maximum height above the stel- 
lar surface of ~ 0.3 stellar radii). The length derived using the 
quasi-static formalism is about 4 times larger, and would thus 
again lead to the "classic" picture of long, tenuous loops. 

Although certainly not small, flaring loops of i ~ 0.5 i?* 
are by no means exceptional, even by the relatively modest 
solar standards. The analy sis of the large flare on Algol seen 
by the SAX observatory ( Schmitt & Favata 1999 ; Favata & 
Schmitt 1999) shows that the picture of large tenuous loops 
implied by the results of the quasi-static analysis can be quite 
misleading, and that, at least in that case, the geometrical size 
of the flaring plasma constrained by the light-curve eclipse 
is significantly smaller than the loop heights derived with the 
quasi-static method. An analysis based on the method of Reale 
et al. (1997) yields, also in that case, loop dimensions which are 
significantly smaller than the ones implied by the quasi-static 
analysis. The much more compact loop size derived through the 



Reale et al. (1997) method is linked with the presence of sig- 



nificant plasma heating during the flare decay, as the implied 
small, dense loop has a very short thermodynamic decay time 
(rth — 400 s). No large diffuse plasma structures at several 
stellar radii from the surface are needed to model the flaring 
region, and a more appropriate picture appears to be one of a 
rather compact, high-pressure plasma structure, whose decay 
is completely dominated by the time-evolution of the heating 
mechanism. 

The present large EV Lac flare shows several characteris- 
tics in common with other large, well observed stellar flares 
discussed in the literature. The light curve has a "double ex- 
ponential" decay, with the initial time scale being more rapid, 
and a slower decay setting in afterwards. A very similar de- 
cay pattern is observed in the large Algo l flare seen by SAX 
as well as on many large solar flares (see Feldman et al. 1995 



for an example). The best-fit metal abundance for the flaring 
plasma also shows what by now appears to be a characteristic 
behavior, i.e. it increases in the early phases of the flare, peak- 
ing more or less with the peak of the light-curve, and then it 
decreases again to the pre-flare levels. In the case of the SAX 
Algol flare the long duration of the event and the high intensity 
of the flare make it possible to show that the metallicity decays 
to its pre-flare value on faster time scales than either the plasma 
temperature or the emission measure. The coarser time cover- 
age of the present flare and the shorter duration of the event do 
not allow such detailed assessment. 

The heating mechanism of the solar (and stellar) corona re- 
mains in many respects an unsolved puzzle, and even more the 
mechanism responsible for large flares. However, it is by now 
rather clear that most sizable flares cannot be explained with a 
simple picture of a (mostly) impulsive heating event followed 
by decay dominated by the free thermodynamic cooling of the 
plasma structure. On the contrary, the evidence from the re- 
cent flow of well studied flare data (including the one in the 
present paper) is that the decay of large flares is dominated by 
the time evolution of the heating mechanism. Thus, the double 
exponential decay observed here as well as in other large solar 
and stellar flare is likely to be a characteristic of the heating 
mechanism rather than one of the flare decay. 

The interpretation of the increase in best-fit plasma metal- 
licity during the flare's peak is still unclear. If the presence of 
a fractionation mechanism is accepted, which causes differ- 
ences in the metal abundances in the photosphere and in the 
corona, the abundance increase during the flare can plausibly 
be explained as due to the evaporation of photospheric plasma 
during the early phases of the flare, on time scales faster than 
the ones on which the fractionation mechanism operates. Since 
the coronal plasma shortly after the impulsive heating is al- 
most entirely made of material evaporated from the chromo- 
sphere, if this scenario were correct the observations presented 
here would imply that the chromospheric metallicity should be 
about three times the coronal one in quiescent conditions. 

Recent calculations (G. Peres, private communication) 
show that the plasma in a flaring loop such as the one respon- 
sible for the EV Lac flare discussed here is not optically thin 
for the strongest lines. This is in particular true for the Fe XXV 
complex at ~ 6.7 keV, which drives the determination of the 
metallicity of the flaring plasma. However, optical thickness 
effects would in general depress the strong line, leading to a 
lower metallicity estimate, and cannot therefore explain the 
metallicity increase observed during the flare. Another possi- 
ble bias to the best-fit metallicity can derive from the thermal 
structure of the flaring plasma, which is not isothermal, even if 
it's being fit with an isothermal model. To assess the magnitude 
of this effect we analyzed the synthetic spectra produced with 
an ad hoc hydrodynamic simulation of a flaring loop, showing 
that the effect is small (< 30%) in comparison with the ob- 
served magnitude of the change (a factor of ~ 3). 

If the heating mechanism responsible for the present flare 
is essentially due to some form of dissipation of magnetic en- 
ergy, an obvious question to ask is what field strength would 
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be required to accumulate the emitted energy, and to keep the 
plasma confined in a stable magnetic loop configuration. A re- 
lated question is whether the present flare could be interpreted 
with a reasonable scaling of the conditions usually observed 
in the solar corona, or whether a different configuration and/or 
mechanism for energy release need to be invoked. Our flare 
analysis allows to make some relevant estimates, under the as- 
sumptions that the energy release is indeed of magnetic origin 
and it occurs entirely within a single coronal loop structure, 
with the characteristics inferred from the analysis of the flare 
decay. An estimate of the minimum magnetic field B necessary 
to produce the event can then be obtained from the relation: 



- even smaller than dictated by equipartition arguments - al- 
ter the convective instability conditions (Ventura et al. 1998) 



X Mo 



(11) 



where _Etot — 7 x 10 erg is the total energy released 
(Sect. 3.3), Bo is the magnetic field surviving the flare and 



Vioop is the total volume of the flaring plasma. As an estimate of 
i?o we take the magnetic field necessary to maintain the plasma 
confined in a rigid loop structure along the whole flare, thus im- 
plicitly assuming that the loop geometry does not change dur- 
ing the flare. From a plasma density n ~ 2 x 10^^ cm~^, and 
a temperature T ~ 100 MK the estimated maximum plasma 
pressure is 6 x IC* dyn cm^; in order to support such a pres- 
sure, a field of ~ 1.2 kG is required. Hence, the total minimum 
magnetic field required to explain the flare is, from Eq. [ll], 
B ~ 3.7 kG, a value compatible with the average magnetic 
field of 3.8 kG, with a surface filling factor of about 60% (and 
evidence for field components of up to ~ 9 kG), measured on 



EV Lac by Johns-Krull & Valenti (1996) at photospheric level 



Although we have used the loop volume in the derivation 
of B, this is not to say that the field fills up the whole volume. 
Rather, our estimates can be interpreted in the framework of 
the flare energy stored in a magnetic field configuration (e.g. a 
large group of spots) with a field strength of several kG, cover- 
ing a volume comparable to the one of the flaring loop. What 
rests as a matter of speculation is how often such a large energy 
release may occur, or in other words, what are the conditions 
required to accumulate such large amounts of magnetic energy, 
especially when the photosphere is so permeated of magnetic 
fields as shown by Fohns-Krull & Valenti (1996). 

In any case, the above scenario suggests the presence of 
large-scale, organized magnetic fields. This is somewhat in 
contrast with the hypothesis that EV Lac is a fully convective 
star, whose activity is powered by a turbulent dynamo, which 
would be expected to produce small-scale magnetic fields. 
Most dynamo theories suggest (Durney et al. 1993) that less 



magnetic flux should be generated by a turbulent dynamo (as 
compared to the case of the solar-type "shell" dynamo) because 
there is no stable overshoot layer where the fields can be stored 
and amplified, and only small-scale magnetic regions should 
emerge uniformly to the surface, because the crucial ingredi- 
ent is small-scale turbulent flow field, rather than large-scale 
rotational shear. 

On the other hand the presence of a magnetic field may sub- 
stantially modify the stellar interior structure. Magnetic fields 



and thus likely modify the structure of the convective enve- 
lope. At the same time convection has the tendency to pump 
magnetic fields downward ("turbulent pumping", Brandenburg 
et al. 1992; Tobias et al. 1998), so that - in a fully-convective 
star - fields may accumulate near the center. Hence, magnetic 
fields are likely to be an important (and thus far essentially un- 
accounted) term in determining the actual stellar structure, and 
any realistic calculation at the low-mass end should consider 
the dynamo-generated magnetic fields as an essential part. In- 
deed, the possibility that a strong magnetic field can lead to the 
formation of a radiative core has been discussed by Cox et al. 
(1981), and this may be the seed for a resurrection of a "shell" 
dynamo mechanism. 
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Appendix A: Physical cliaracteristics of EV Lac 

EV Lac (Gl 873) is a dMe dwarf, classified as M3.5 (Reid 
et al. 1995), at a distance (from the Hipparcos-measured par- 
allax) d — 5.05 pc. It is considered a single star, with no 
evidence of companions, and is a slow rotator, with a pho- 
tometrically determined rotation period of 4.378 d (Pettersen 
et al. 1992). The projected equatorial velocity has been de- 
termined from Doppler broadening of the spectral lines at 
vsini = 4.5 ± 0.5 km s"^ (pohns-Krull & Valenti 1996| ) and 
V sini = 6.9 ± 0.8 km s~^ ( pelfosse et al. 1998 ); this rotation 



V sini < 5 km s~^, and a tail of rare fast rotators with veloci- 
ties of up to ~ 50 km s^^. EV Lac lies at the upper end of the 
slow-rotator distribution. 

The absolute magnitude is My = 11.73, which, given a 
color index {R — I)c = 1-52 translates in A/boi = 9.40 
( Delfosse et al. 1998 ), corresponding to Lboi ~ 5.25 x 
10'^^ erg s~^. Using the mass-luminosity relationship of 



Baraffe et al. (1998) and the /f-band luminosity (to be pre- 
ferred given the independence of the relationship between 
mass and j^T-band lum inosity from metallicity) Mk = 6.78 
(Delfosse et al. 1998) we derive a mass of ~ 0.35 Mq. No 



velocity can be reconciled with the observed starspot modu- 
lation period only if the inclinat ion is high (> 60 d eg). The 
rotational velocity of M dwarfs ( pelfosse et al. 1998 ) is char- 
acterized by the bulk of them having a narrow distribution with 



photospheric abundance analyses of EV Lac are known to us, 
although Fleming et al. (1995) report a near-solar metallic- 
ity based on broad-band photometry. The radius for a solar- 
metallicity 0.35 Mq dwarf is, from the models of Chabrier & 
Baraffe (1997), i? ~ 2.5 x 10^" cm, or i? ~ 0.36 Rq, (assum- 
ing the star is old enough, given that such a low-mass star will 
contract gravitationally until it's ~ 300 Myr old). 

Stellar structure models show that the radiative core of low- 
mass stars shrinks with decreasing mass, disappearing com- 
pletely in mid-M dwarfs, so that late-M dwarfs are expected 
to be completely convective. Solar-type dynamos are thought 
to require an interface between the convective envelope and the 
radiative core (the a-il shell dynamo model) and are thus not 
expected to be present in the cooler M dwarfs. Given that how- 
ever X-ray activity is present and common even in pu rportedly 
fully convective late M dwarfs ( Barbera et al. 1993 ; Schmitt 
et al. 1995) a different type of dynamo mechanism must be 
operating in them; it has recently been suggested that small- 
scale magnetic fields can be generated in convection zones by 
a turbulently driven dynamo ( purneyetal. 1993 ). This - which 
could be also at work in higher mass stars with varying degrees 
of efficiency, depending on the rotation rate - would therefore 
provide the only magnetic field generation mechanism in fully 
convective low-mass dwarfs. 

The predicted mass at which stars become fully convec- 
tive dep ends on the details of the ph ysics adopted in the stellar 
models. Chabrier & Baraffe (1997) use non-grey atmospheres 
to put the fully convective limit at 0.35 Mq (i.e. just at the 
estimated mass of EV Lac), independent of metallicity in the 
range 0.01 x Zq < Z < Zq. Different calibrations for the 
mass-luminosity relationships (see discussion in Delfosse et al. 
1998) push the fully-convective limit toward somewhat higher 
masses, so that EV Lac is most likely fully convective, and thus 
its activity is likely to be driven purely by a turbulent dynamo, 
which is not, among other things, expected to have a solar- 
like cyclic behavior, and which is expected to generate mag- 
netic fields on a smaller spatial scale (related to the scale of 
the turbulent flow fields) than an a-ft dynamo. The picture is 
complicated by the fact that strong magnetic fields may influ- 
ence the stellar interior structure maintaining a radiative core 
at masses lower than the theoretical limits for spherically sym- 



metric, non-magnetic stars (Cox et al. 1981), and thus the mag 



netic fields of very active, low-mass stars could still be (par- 
tially) powered by an a-il shell dynamo. 
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A.l. Previous X-ray observations 

The high activity of EV Lac had been noted well before the 
advent of high-energy observations from its large optical and 
UV flar ing rate, with some truly excepti onal optical flares ob- 
served: Roizman & Shevchenko (1982) report the occurrence 
of a 6.4 mag J7-band flare, lasting 6.4 hr, with a peak luminosity 



of ^ 



10^2 erg s-i 



and a total energy output of > lO'^^ erg. 



EV Lac was first observed as a 



quiescent soft X-ray 
source by EXOSAT ( [Schmitt & Rosso 1985| ) - although it had 
been detected earlier by HEAO-1 during a strong flare with 
logix = 28.7 erg s^^. It was not observed by the Einstein 
observatory, while it was observed by ROSAT both in the All 
Sky Survey (RASS - in which a major flare was also detected) 
and in pointed mode. The R ASS quiescent X-r ay luminosity 
was logix = 29.08 erg s^^ ( Schmitt et al. 1995 ), correspond- 
ing to logLx/^boi — —2.6. It was also the subject of one 
SAX and several PSPC pointings, analyzed in Sciortino et al. 
(1999), during which its coronal emission shows both a contin- 
uous variability of about a factor 2-3 and the occurrence of an 
intense flare, with an increase of the emission in the PSPC of 
about 10-fold. 

The HEAO A-1 Sky Survey experiment (Ambruster et al. 
1984) detected two flares from EV Lac, with peak X-ray lumi- 
nosities (in the 0.5-20 keV band) of logLx ~ 29.5 erg s^^ 
and logLx = 30.3 erg s^^. The most energetic of the two 
flares represented a peak enhancement of ~ 50 over the qui- 
escent X-ray luminosity (Lx — lO^^'^ erg s^^), and its de- 
cay e-folding time was roughly estimated to be of the order of 
< 5 ks (although the very sparse temporal coverage prevents an 
accurate determination of the flare's decay). Ambruster et al. 
(1984) have also estimated the physical parameters of the two 
flares (assuming a peak temperature of ~ 2 x 10^ K), resulting, 
for the smaller of the two flares, in a peak emission measure 
EM ~ 2 X lO^'^ cm^'^, a density n ~ 5 x 10^^ cm^'^ and 
a loop length L ~ 5 x 10^ cm. The magnetic field strength 
necessary to confine the plasma was estimated at 200 G. Only 
the peak emission measure is reported for the second flare 
(EM ~ 2 X lO^'^ cm^'^). A small flare was seen in one of 
the EXO SAT observations, in the LE detector, as discussed in 
detail by Pallavicini et al. (1990) . Its rise time (1/e time) was 
~ 600 s, and its decay time was ~ 4.5 ks. At peak the flare rep- 
resented an enhancement of only ~ 3 times over the quiescent 



X-ray flux, with a peak X-ray luminosity Lx 
and a total energy released in the X-rays E 2 



10^ erg s ^ 
10^2 erg. The 

lack of energy resolution of the EXOSAT LE detector made it 
impossible to perform an analysis of the flare's decay. 

Schmitt (1994) derived loop parameters for the long RASS 
flare by fitting the flare decay parameters within the framework 
of the quasi-static formalism of van den Oord & Mewe (1989). 
The maximum observed temperature is T ~ 30 MK, the de- 
cay timescale is r ~ 38 ks, and the peak emission measure is 
EM ~ 1.5 X 10^2 cm~'^. The loop length derived through a 
quasi-static analysis is large, at L ~ 6 x 10^^ cm ~ 10 i?, (with 
an inferred flaring volume 1/ ~ 3 x 10^^ cm'^) and the plasma 
density is correspondingly low, at n ~ 3 x 10^° cm^^. The to- 



tal thermal energy was estimated at ~ 9 x 10'^^ erg. Schmitt 
(1994) also analyzed the EV Lac PSPC flare decay using the 



two-ribbon model of Kopp & Poletto (1984) 



